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Abstract 

We present spline wavelets of class C n (R) supported by sequences 
of aperiodic discretizations of R. The construction is based on mul- 
tiresolution analysis recently elaborated by G. Bernuau. At a given 
scale, we consider discretizations that are sets of left-hand ends of 
tiles in a self-similar tiling of the real line with finite local complex- 
ity. Corresponding tilings are determined by two-letter Sturmian sub- 
stitution sequences. We illustrate the construction with examples 
having quadratic Pisot-Vijayaraghavan units (like r = (1 + \/5)/2 
or t 2 = (3 + v5)/2) as scaling factor. In particular, we present a 
comprehensive analysis of the Fibonacci chain and give the analytic 
form of related scaling functions and wavelets. We also give some 
hints for the construction of multidimensional spline wavelets based 
on stone- inflation tilings in arbitrary dimension. 
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1 Introduction 



Under the name "wavelet" is commonly understood a function if)(x) G L 2 (IR) 
such that the family of functions ipj,k(x) = 2 j / 2 ip(2 j x — k) for j,k G Z forms 
an orthonormal (or at least a Riesz) basis for L 2 (R). A function generat- 
ing through dilatations and translations an orthonormal basis for L 2 (IR) can 
be found through a "multiresolution analysis of L 2 (IR)" (shortly MRA), a 
method settled by S. Mallat [1]. The dilatation factor is usually 9 = 2. In- 
deed, the construction of a wavelet basis within the MRA framework relies on 
the fact that the lattices 2~ j Z are increasing for the inclusion. This property 
is preserved only when 9 is an integer. Then one can raise a natural question: 
what about choosing another number 9 as a scaling factor? A first answer is 
given in the work of P. Auscher [2] . The following problem was considered: 
given a real number 9 > 1, does there exist a finite set {ipi,ip2, ■ ■ ■ ,ipe} of 
functions in L 2 (R) such that the family 9 j/2 ipi(9 j x -k), j,keZ, 1 < i < £, 
is an orthonormal basis for L 2 (IR)? The author then proved that a basis of 
this type exists if 9 is a rational number. More precisely, for 9 = p/q > 1, p 
and q being relatively prime integers, there exists a set of p — q wavelet func- 
tions satisfying the previous condition. It was still an open question whether 
there exists another generalization of wavelet basis with an irrational scaling 
factor. On the other hand, in 1992 Buhmann and Micchelli [3] proposed 
a construction of a wavelet spline basis corresponding to non-uniform and 
non-self-similar knot sequences. Further [4, 5, 6] studies have been recently 
devoted to this problem in higher dimension, like answering the question of 
characterizing functions ip, dilation sets V, and translation sets T, such that 

{\det(D)\^(Dx - X)\D G £>, A gt| forms an orthonormal basis for L 2 (R) 

[6]. As a matter of fact, it was proved in [5] that for any real expansive (i.e. 
all eigenvalues have modulus greater than one) d x d matrix A, there exists a 
measurable E C R d such that {\det A\ j/2 ^ E (A j x - I), j G Z, / G Z d ] forms 
an orthonormal basis for L 2 (M. d ). 

In 1996, a Haar wavelet basis with an algebraic irrational scaling factor 
and only one generating wavelet ip(x) was given in [7]. This wavelet ba- 
sis of L 2 (R) lives on the nested sets t ? Z t having the following structure: 
T'jl 2 i\){jix — A), j G Z, A G A, where r = (1 + v / 5)/2 is the golden mean. Set 
Z T means set of "tau-integers" and set A is the set of admissible translations. 
One outcome of the latter work is that admissible translations are generically 
irrational for an irrational scaling factor. This seems to be a common feature 
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to aperiodic translational sets. In 1998 G. Bernuau [8, 9] settled a general 
construction of spline wavelets living on locally finite (more precisely with 
finite local complexity) and self-similar Delaunay (or Delone) sets. His ap- 
proach was mainly inspired by important results previously obtained by P.G. 
Lemarie-Rieusset [10] in the more general framework of stratified nilpotent 
Lie groups. We learn from these last works that irrational factors combined 
with specific properties of sequence of discretizations imply finitely many 
generating wavelets and that we also need more functions named "scaling 
functions" in an appropriate modification of MRA. Actually the number of 
spline scaling functions and spline wavelets does not depend only on the 
scaling factor but also on the polynomial order of these functions. 

In a previous letter [11], we have presented a definition of multiresolution 
analysis for an infinite sequence ••• C jF/r 2 ^ 2 C J-'/t 2 ^ C J-/t^ +2 C ••• 
of aperiodic discretizations of R. Corresponding wavelets have been defined 
and the elementary Haar example was given as an illustration of the method. 
Starting from the Haar case, one could be attempted to explore less trivial ex- 
amples in which compact support and orthogonality are required for wavelet 
families of class C n , n > 0. This is a strong constraint for such Daubechies- 
like wavelets which live on aperiodic discretizations, and their existence is 
not guaranteed at the moment. Our present strategy is to drop out the or- 
thogonality condition and to rather explore certain well-known quasiperiodic 
counterparts of the dyadic spline wavelets. In the present paper, we shall 
closely follow the procedure rigorously settled by G. Bernuau. Let us briefly 
summarize the general setting of the Bernuau construction. Let A C Rbe 
a Delaunay point set in the real line. By Delaunay we mean that A is uni- 
formly discrete (the distances between any pair of points in A are greater 
than a fixed r > 0) and relatively dense (there exists R > such that K. is 
covered by intervals of length 2R centered at points of A). In addition to 
this Delaunay structure, we demand that the set A be 

1. self-similar: there exists a number 9 > 1 (inflation factor) such that 



6A C A, 



(1) 



2. with finite local complexity: for all R > 0, the point set 




(2) 
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is finite. This means that local environments of points in A are not 
different in infinite fashions. 



Typically, such sets A are mathematical models for one-dimensional struc- 
tures having a long-range order, like quasicrystals. Our aim here is to con- 
struct a Riesz basis of L 2 (R), the elements of which are of the affine wavelet 
type: 

x — k), j G Z, k G 9 1 A, succ(k) ^ A, (3) 

where succ(k) is the nearest right neighbour of k in the set # _1 A, and for 
which the set {if) K (x)} of mother wavelets is finite and the simplest possible. 

At this point, we recall that (v n ) is a Riesz basis of a separable Hilbert 
space V if and only if each v G V can be expressed uniquely as v = Y2 n a nV n 
and there exist positive constants A and B, < A < B, such that 

II 2 

A^2\a n \ 2 < ||^a n v n <Bj2\ a n\ 2 (4) 

ran n 

for all sequence of scalars a n . We can say that the i> n 's are strongly linearly 
independent. 

In the next section, we shall present a survey of preliminary results con- 
cerning the space of splines based on A and at the heart of the construction 
of wavelets and corresponding multiresolution analysis. In Section 3, we re- 
call the Bernuau theorems about the existence and the characterization of 
the wavelet basis itself. Section 4 is devoted to the description of Delau- 
nay sets based on two-letter Sturmian substitution sequences. These sets 
have as scaling factor special algebraic integers, with generic symbol f3, and 
named quadratic Pisot-Vijayaraghavan units or more simply UPV2. Among 
the latter one finds those numbers which are of interest in quasicrystalline 
studies: r = (1 + \/5)/2, r 2 = (3 + V / 5)/2 (for pentagonal and decagonal 
cases), 00 = 1 + \/2, oo 2 = 3 + 2v^2 (for octogonal case), and 5 = 2 + 
(for dodecagonal case). In this context, we shall consider two types of point 
sets: the /5-integers and certain model sets included in the latter. Original 
results are given in Sections 5-7. Haar wavelets living on /9-integer sets with 
corresponding scaling equations are given in Section 5. Section 6 is devoted 
to the study of one important example of model set, namely the Fibonacci 
chain, from a lexicographical point of view (see Proposition 6 and Properties 
7-9) and prepares the section 7 in which related splines and wavelets together 
with their scaling equations are explicitly constructed. In Section 8, we con- 
sider a class of self-similar tilings of M. d having the so-called stone-inflation 
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property. For instance, we can find stone-inflation tilings among Penrose or 
triangle tilings of the plane. We just sketch the method of construction of 
Bernuau spline wavelets adapted to such tilings. 



2 Space of splines for multiresolution analysis 

2.1 The definition of the space Vq (A) 

Any Delaunay set A determines a space of splines of order s, s > 2, in the 
following way 

Definition 1. Let s > 2. Then Vq^(A) is the closed subspace of L? 
defined by V^{K) = G L 2 (R) | ^f(x) = £ AeA a A <f, 

An equivalent definition is given in terms of the restriction of func- 
tions to intervals determined by consecutive elements of A. Suppose the 
latter is defined by the increasing one-to-one map from Z : Z 3 n — > 
A n , • • • A n _i < A n < A n+ i • • • . There results from Def. 1 that Vq^(A) = 

j/ G C s ~ 2 , f G L 2 (M,), f [Xn x +i] is a polynomial of degree < s — 1 j. There- 
fore, Vq S ^(A) is the space of splines of order s with nodes in A. Let us now 
give a classical result about the existence of a Riesz basis for V^(A) [12]. 

Theorem 1. For all Delaunay sets A = {X n }nei C K. and for all s > 2, 

there exists a Riesz basis {B^\\ G A} o/Vq S ^(A). The function (called 
B-spline) is the unique function in Vq^(A) satisfying the following conditions: 

(i) supp B ( x s) = [A, A'], where A' G A. 

(ii) The interval (A, A') contains exactly s — 1 points of A. 
A'- A 



(Hi) j B^ = 
Jr 



See [12] for proof. Note that (i) and (ii) give precise information on the 
(compact) support of whilst (iii) is a normalization condition. 
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2.2 Construction of B-splines of order s, B^\x) 



The construction of B^'(x) can be carried out in various ways: 

• by recurrence [13], with the use of the formula valid for all A n G A 



1 if A n < x < A n+ i, 
otherwise, 



(5) 



B%(x) = u^B^ix) + (1 - cu s , n+ i(x)) Bl'Wx), for s > 2 (6) 
with u 8 , n (x) = {x- \ n )/ (A„ +S _i - A n ), 

by using the condition of minimal support, which means that V/ G 
(A), 

supp / C supp B ( x s) =>• / oc Bj^, 



by inverse Fourier transform which gives the explicit form 



B^ (x) oc Fourier 1 



(7) 



(8) 



in which {ai(n)} is the unique solution of the linear system 

s 

J2(*n+l) j Xl = 0, < j < S ~ 1, 

(-!)• 



1=0 



(9) 



■n+zj #z 



z=o 



Note that the above linear system is easily solved since it involves Vander- 
monde determinants: 



ai(n) 



-1) 



s\ 



j~J (A n +Z — A n+ //) I J (A n+ ;/ — X n +l) 



.0<l'<l 



Kl'<s 



(10) 



From (8) one can check that the s th derivative of B^\x) is given as a finite 
linear superposition of Dirac masses located at the points \ n +i, < I < s: 

s 



dx 
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Since the Fourier transform of B^'(x + X n ) depends on the s-tuple (A n+ i — 
A n , . . . , A n+S — A n ) only and since such s-tuples assume their values in a finite 
set for n varying in Z if A has finite local complexity, we can assert the 
following: 

Proposition 1. Let A C H. be a Delaunay set of finite local complexity. Then 
the B-splines of order s based on A are of the form B^\x) = <f)\(x—X), A G A, 
where the set {(f)\(x), A G A} is a finite set of functions with compact support. 

Therefore, in the finite local complexity case, it is possible to partition 

q 

the indexing set Z for A into a finite set of equivalence classes n, Z = U_n 

n=0 

where class n is defined by 

n = {fc G Z | B[ s) k (x + X k ) = B { £(x + X n ) Vx} . 

q 

Correspondingly, for a given s, the point set A is partitioned into A = O As 

n=0 

with An = {Xk G A | k G n}. The relation A^ G A n means that A^ and X n are 
left-hand ends of identical s- letter words if we identify each interval (Xk, Xk+i) 
with a letter of the allowed alphabet. To each class n is biunivocally asso- 
ciated the function (f) n (x) = 4>\ k (x) = B^(x + Xk), k G n. In this way, the 

space Vq S \A) decomposes into the direct sum 

^o W (A) = 0^, (11) 

n=0 

where V 0jn is the closure of the linear span of the functions <p n (x — X k ), k G n. 

2.3 Self-similarity and mult iresolut ion analysis 

Let A be a Delaunay set of finite local complexity and self-similar with infla- 
tion factor 9 > 1 (9 A C A). Changing the scale allows us to define subspaces 
Vf } (A), j G Z. 

Definition 2. vj s) (A) = {f(x) G L\R) \ £,f(x) = E A£A ^a} ■ 

We now have at our disposal an inductive chain of spaces allowing analysis 
at any scale. More precisely, with the above notations, 
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Proposition 2. The sequence of sub spaces (Vj (A))j e z is a -multiresolution 
analysis of L 2 (Jit), i.e. 

(i) for any j is a closed subspace of L (M), 

(ii) ■■■ CV^(A) cV {s \A) cV{ s \A) c ■■■ , 
(™) U je z V j S \ A ) ls dense m l2 ( r ); 

m n i6Z ^ w (A) = {o}, 

(v) f(x) G vj s) (A) if and only iff(9^x) G V^(A), 

(vi) there exists a finite number of functions <pn(x) G Vq S \ called scaling 
functions such that {<f)n(x — ^k)}k<=n o<n<q ^ s a R^ esz basis in Vq S \ 

The proof is straightforward from Definitions and Proposition 1. 

3 Spline wavelet basis 

As is well known, the whole wavelet basis is obtained by scaling the elements 
of a wavelet subfamily living at a given scale of the multiresolution. For 
convenience, we choose here the O th -scale and this subfamily is precisely a 
basis of the orthogonal complement W^(A) of V^? (A) in Vq (A): 

V (s) (A) = Vi? (A) © ± W { _ S }(A). (12) 

This decomposition can be interpreted as the building up of the "(j = 0) th - 
scale" content of the multiresolution (i.e. Vq^) from the adding of necessary 
"details" (i.e. W^l) to the existing content at the next larger scale (i.e. 
V^i). A crucial step in the characterization of these "details" is the following 
existence theorem. 

Theorem 2. (Bernuau [8, 9] Let A e R be a Delaunay set of finite local 
complexity, self-similar with scaling factor 9 > 1. Then, for all s > 2, there 
exists a Riesz basis of W^l (A) of the form 

{([% ~ An), A n G A, A n+1 i 9A] , where 
is a finite set of functions with compact support. 
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We here sketch the proof by just listing a sequence of intermediate results 
given in [8, 9]. A preliminary characterization of W*f\ is given by the following 

proposition. Let us recall that H S (M) denotes the Sobolev space of functions 

d a 

f G L 2 (R) such that -j-^f, with a < s, is element of L 2 (R). 

Lit/' 

Proposition 3. For s > 2, W^l is the space of s th derivatives of elements 
in H 2s (M.) which vanish on OA: 

W^}(A) = ^—K 2s = I ^—h. h e H 2s (R) | VA G 9A, //(A) = 

Proof. As indicated in Equation (12) let us determine the orthogonal comple- 

d s 

ment of V^(A) in V (s) (A). Let / G V^(A) and h = —h u where h x G K 2s . 

ax s 

Then 

Adx = ( /,il ;ftl \ = (-ir(il ; /, ftl ) = o. 



/ 

Jr 



Conversely, let / G L 2 (1R) and orthogonal to W*fl (A) = — —K 2s . Let us 

\XJb 

choose a function (p G C£° such that <p(0) = 1 and with such small compact 
support that its translates ip(x — A), A G OA, have disjoint supports. For all 
g G Cg°, #(x) - J2xee\9W^( x ~ A) is in fT 2s . Accordingly 



/ 



( ^ " E - A) ) dx = 0, 

Ae6»A 



so we have 



dx s / \ dx 

Ae6»A 



which means that / G V?{ (A). □ 

Remark 1. ^4n equivalent characterization is to write that, fors > 2, wl^A) 

d s 

i/ie set of functions f G £ 2 (IR), / = -z — h, for which h G Vq 2s ^(A) and 



h\ 0A = 0. 



dx s 
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In view of this remark, the explicit construction of these wavelets neces- 
sitates the introduction of the following space 

w[ s l(A) = {feV {s \A)\f\ eA = 0). 

This is a closed subspace of Vq^(A). If we consider in particular W^i\A), 

then, by derivation of its elements, we will obtain functions in w!f)(A). 

dx s 

Now, we need functions in wi^(A) with minimal support. Let us consider 
the following subset of integers: 

E = {neZ | A n+ i £ OA} . (13) 

Next, for all n G E, let us define the unique function \t4 s ^ G (A) satisfying 
the following conditions: 

(i) supp^n"* is compact and included in [A n ,oo), 

(ii) ¥:\\ n+1 ) = i, 

(iii) / G W[ S }(A), supp/ C supp^l s) / cx *i s) . 
Then we have another important result. 

Theorem 3. The set |^n\ n G -E 1 j zs a i?zesz fraszs o/ £/ie space W 7 "!^ (A). 

The proof of this theorem is also given in [8, 9] and goes through a list of 
intermediate properties which we give here without proof: 

Property 1. For all n G E, the support of is the interval [A n ,A n+ jv n ], 
where N n is the smallest number equal or larger than s such that 
N n = s + #{(\ n ,\ n+Nn ) n 9 A}. 

Note that N n is also defined as the smallest number such that the interval 
(A n , A n+ Ar n ) includes exactly s — 1 points from A \ 6A. 

Property 2. There exists a finite set Q s of functions with compact support 
such that, for all k G E, ^ (x + A&) G Q s . 



Property 3. All functions from W^l(A) with compact support are linear 



combination of the ^u^, k G E. 
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Property 4. The set of functions in W^i(A) that have compact support is 
dense in w[ s }(A). 

Property 5. For all m G Z, define I m as the subset of k G E for which 
is not identically equal to zero on the interval [A m , A m+ i]. Then s — 1 < 
#(/ m ) < s. 

Property 6. There are two constants < A < B , depending on s and A 
only, such that, for any sequence (ak)keE with finite support, we have 



Aj2\a k f<\\J2a k ¥:\<Bj2 

keE k<=E ~ keE 



1 2 

" \a k \ 



Note that if A = 1 = B then the set , n G i?| is an orthonormal 

basis. The main result of this section is a consequence of all the above 
statements: 

Theorem 4. (Bernuau [8, 9]) Let A = {A n } neZ C R be a Delaunay set 
of finite local complexity, self-similar with factor 9 > 1. Let us denote the 
elements of 9~ 1 A by K n ,n G Z, /t n = # _1 A n . T/ien /or a// s > 1 i/iere exisfo 
a Riesz basis of L 2 (R) o/ the form: 

{e j / 2 ^}(e j x - Kn), K n g e-'A, Kn+1 $ a, j g z} , (14) 

where | a /imte set o/ compactly supported functions of order C s ~ 2 . 

Proof. For A n G A, define C\^( x ~ — ~^~^n S \x). Then, according to 
Property 2, jcjj} i s a finite set for A n+i ^ OA. These functions together with 
all their admissible translates C\}{ x ~~ ^n) form a Riesz basis of W^(A). We 
know that the sequence of spaces (v} s \A)) is a multiresolution analysis 
of L 2 (R), so we get the orthogonal sum 

L 2 (R) = ± wf(A) 
where /(a;) G W^ s) (A) G W (s) (A). 
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For convenience we shift by one the scale and choose the L 2 -normalization 
in order to define 

*- w = M e <)(A) ' (15) 

Then we can assert that the set ^6^ 2 ^k}(9^x — n n ), n n G -1 A, /t n+ i ^ a| 
is a Riesz basis of W^ s) (A). By union on all j G Z we obtain the theorem. □ 

4 Self-similar discretizations of R with UPV2 
scaling factor 

4.1 Two-letter substitution sequences, quadratic PV 
numbers, and bet a- integers 

Let a be a positive integer and consider the following two types of two-letter 
substitution (L for "long" and S for "short") 

a times 

q : { L -> L^LS, with a > 1. (16) 
S -> L, 

a— 1 times 

" <, " tiroes ' with ° > 3 - ( 17 ) 

S -> L^LS, 

The corresponding substitution matrices read respectively ( f J ) and ( "Ig i ) • 
The characteristic equation for the former is 



X 2 - aX - 1 = with roots <j Z^ _\- pl = a ^ ( 18 ) 



For the second matrix we have 



= P, 



X 2 -aX + l = with roots <| _ jl'= p' — q /J. (19) 

In both cases, the largest root /3 is > 1 whilst its Galois conjugate j3' G 
(— 1, 1). For this reason, one says that /5 is a quadratic Pisot-Vijayaraghavan 
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and a unit for being invertible in their respective extension ring Z[/3] (the 
lowest-degree coefficient in (18) and (19) is ±1). Henceforth, we shall denote 
by B- (resp. B+) the set of UP"VYs obeying (18) (resp. (19)). 

Let us now associate to the letter L a tile of length 1 and to S a tile 
of length 1/(3 for the case (16) and 1 — 1/(3 for the case (17). Starting 
from the origin of the real line with L on the right, we apply to it the 
substitution The set of nodes of the resulting tiling of the half-line 
is called the set of positive "beta- integers" and is denoted by Zt. Taking 
the symmetric — Z^ of Z^ with respect of the origin, we obtain the set of 
"beta-integers" Z^ = Z^ U (— Z^). This set is clearly a Delaunay set which 
is of finite type complexity and self-similar with factor (3: fiZp C Zp. By 
construction, it is also symmetrical with respect to the origin: — Z^ = Z^. 
This name of beta-integers comes from the fact that if we retain in the set 
of real numbers all those numbers which are polynomial in (3 when written 
in "basis then they form a totally ordered discrete set that precisely 
coincides with Z^. We recall here that the writing of a real positive number 
x in irrational basis (3 > 1 means a unique /5-expansion x = X^=-oo x ifl l = 
XjXj—i ' ' ' X\Xq ' x —\ ' ' ' X\ ' ' ' m which the expansion coefficients assume their 
values in the alphabet {0, 1, . . . , \_(3\ integer part of (3 } in agreement with 
the existence of allowed words determined by the so-called greedy algorithm 
(for more precisions, see for instance [14, 15]). In the "minimal" case a = 1 
of the first category (18), (3 — r = (1 + v / 5)/2 ss 1.618 • - • , the rules are 
particularly simple, since there the alphabet is {0, 1} and the constraint is 
that no two-letter string 11 should appear in any r-expansion words. The 
first pieces around the origin of the associated tiling are shown in the picture 
below. 

S LLSLLSLLS 

e e e e e e e e e 

-r 2 -l -t 2 -t -1 It t 2 t 2 + 1 

Figure 1: Tau-integers around the point 0. 

Note that the case coming just next this minimal one, namely a = 2, is 
the "octogonal" number (3 — uj — 1 + y/2. 

The minimal case in the second category corresponds to a = 3 and yields 
the UPV 2 = (3 + V5)/2 = r 2 = 1 + r w 2.618 - • •. The alphabet is 

n times 

now {0, 1, 2} and the constraint is that no two-letter string 21 • • • 12, n > 0, 
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should appear in any r 2 -expansion words. The first pieces of the associated 
tiling on the right of the origin are shown in the picture below. 

L L LSLLSLSL 

e e e e e e e e e 

1 2 (3 (3+1 (3 + 2 2/3 2/3 + 1 (3 2 

Figure 2: Beta-integers around the point for (3 = r 2 . 



4.2 Model set discretizations of M 

Another way of obtaining self-similar Delaunay sets with finite local complex- 
ity is to resort to the so-called Cut-and-Project method which has become 
like a paradigm in quasicrystalline studies. We shall adopt here the formalism 
set up by Meyer [16, 17]. A 1+1-cut and project scheme is the following 

R RxR R 

U (20) 
D 

where D is a lattice. The projection hi\d is 1-to-l, and 7r 2 (-D) is dense in R. 
Let M = tti(D) and set * = 7r 2 o (tyiId)^ 1 

* : M — ► R. (21) 

The set A C R is a model set if there exist a cut and project scheme and 
a relatively compact set Q C R of non-empty interior such that 

A = {x G M | x* G 0} e Tj n . (22) 

The set Q is called a window. 

As an illustration, we describe one type of Fibonacci chain with scaling 
factor j3 = t 2 . Consider the cut and project scheme (20) with D = Z 2 and 
7Ti(Z 2 ) ~ Z[t] = {a + br \ a, b G Z}. The map (21) is identical up to a factor 
to the Galois ring automorphism x = a + br — >■ x' = a — K An example of 
Fibonacci chain T T [18] is given by choosing the semi-open interval [0, 1) as 
a window. 

/ h 

J r T = {x = a + br\x = a G f2 = [0, 1)} 

T 

= {..., -r 3 ,-r,0,r 2 ,r 3 + l,r 4 ,...}. 
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It is the set of left endpoints of a quasiperiodic tiling of R with 2 tiles L et 
S, of respective length r 2 and r, generated by the substitution rules 

Starting from S to the left and from L to the right we get a biinfinite word 

• • • LLSLS | LLSLLSLS • • • . 

Note that r 2 jF r c T T and that this tiling is, by construction, stone inflation, 
which means that all the tiles when scaled by the factor r 2 can be packed 
face-to-face from the original ones. Therefore, the UPV2 r 2 G B + is the 
scaling factor of the substitution q. The first pieces of this Fibonacci tiling 
around the origin are shown in the picture below. 

L SLSL LSL L S 

© e e e e 9 9 e e 

-r 3 -r -r 3 -t r 2 r 3 + l r 4 r 4 + r 2 r 5 + 1 

Figure 3: Fibonacci chain around the point 0. 

More generally, for (3 G B we shall consider the model set: 

= {x = Z[(3] I x' G Q = [0, 1)} (24) 

where the prime " ' " designates the ring automorphism x = a + b(3 — > x' = 
a + bp'. Let us give a set of properties concerning T$ when (3 G B. Proofs 
can be found in [19, 20, 21]. 

Proposition 4. Suppose that f3 G B-. Then the model set is a self- 
similar Delaunay set of finite local complexity, with scaling factor f3 2 . It can 
be characterized either as a subset of the beta-integers or as set of nodes of the 
stone inflation tiling associated to a substitution sequence. More precisely: 

(i) is obtained from through the sieving procedure 

= {x G Zp I x' G [0,1)}. (25) 

(ii) is the set of left endpoints of a quasiperiodic tiling o/R with 2 tiles 
L et S, of respective length (3 + 1 and (3, generated by the substitution 
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rules 



a times 



a—1 times a—1 times 



a—1 times a times 



L -> LS-SLS-SL- S---S LS---S, 



a—1 times 



(26) 



a— 1 times a—1 times 



a—1 times a times 



[S -> LS--'SLS-SL--- S--SLS--S, 



u^/i a > 1. T/ie taZmg obtained by starting from S-origin-L in both 
directions. 

Proposition 5. Suppose that (3 E B + . Then the model set is a self- 
similar Delaunay set of finite local complexity, with scaling factor (3. It can 
be characterized either as a subset of decorated beta-integers or as set of 
nodes of the stone inflation tiling associated to a substitution sequence. More 
precisely: 

(i) is obtained from the decorated beta-integers TLp = + {0, ±-|} 
through the sieving procedure 

Tp = {xE% | x'e [0,1)}. (27) 



(ii) Alternatively, reads as 



^ u ((-zj) \ {0} + i 



(28) 



fiiij JF^ £/ie set of left endpoints of the quasiperiodic tiling o/R to£/j 2 ti/es 
L e£ S, o/ respective length f3 and {3 — 1, generated by the substitution 
rules (17). The tiling is obtained by starting from S-origin-L in both 
directions. 

Note that the minimal case (3 = r is exceptional in the sense that there is 
equality between both extension rings Z[r] = Z[r 2 ] and so one can deduces 



from (28): T T = (3 



Z+ U ((-ZJ) \ {0} + 



, with /3 



r 2 . 
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5 Haar wavelets for /3-integers 



We now turn our attention on the construction of multiresolution analysis of 
the Haar type based on non- negative /3-integers Z^, the extension to the full 
Z^ being carried out by simple symmetry with respect to the origin. Let us 
denote by space the closure in L 2 (R + ) of the linear span of all positive 
admissible translates of normalized characteristic functions </>l(^) and <fts (x) 
supported by intervals of lengths |L| and |S| respectively More precisely, 

V + = span{0 L (x- A L ),0 s (a:- A s )} ALeA + iAseA +, 

where is the set of left-hand ends of tiles L in Z^ and Ag is the set of 
left-hand ends of tile S in Z^. 

Two possibilities exist for constructing basis made up with piecewise con- 
stant functions. Either one builds orthonormal basis, at the price of increas- 
ing the number of mother functions, or one just requires a Riesz basis and 
then we need two mother functions only. 



5.1 Orthogonal basis of Haar type 

Let us first discuss the case (16). Due to the two-letter substitution rules, two 
possible distances between points exist in the associated tiling, namely |L| = 
1 and |S| = 1/(3. Consistently, two scaling functions exist, one per type of 
tile, and their refinement equations are precisely based on these substitution 
rules. For a > 1 the two scaling functions read: 

<f> L (x) = l m (x), (29) 
Mx-a) =p 1/2 l [a , a+1/f3) (x). (30) 

Consequently, the scaling (or refinement) equations for a > 2 are the 
following: 

a-l 

Mx) = ^ x - l ) + P~ 1/2 MPx - a), (31) 

1=0 

fa(x-a) = [3 1/2 <f) L (fa-a[3). (32) 

The corresponding orthonormal wavelet basis is built from a mother 
wavelets {ip-L,i}iZo ■ One method of construction of these functions is to 
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Gram- Schmidt orthogonalize and normalize in the sense of L 2 -norm the set 
of functions {^(x), 4>l(/3% — 1), • • • , 4>\Xfi x — (a — 1))} living on the 

interval [0, 1]. An orthonormal basis for = V + (B± W + is then obtained 
by collecting together all translates of our orthonormal set at points of 
and all translates of the 0s 's at points of Ag. After removing functions in 
V + we get the basis of W + . 

The case a — 1, (3 — r is particularly interesting. In this case, A^ = i~Z+ 
and A+ = r 2 Z+ + 1. There are two scaling functions but only one wavelet, 
and so the refinement equations read as: 

<Ph{x) = <Pl(tx) + T~ 1/2 (f) S (TX - 1), 

(fis(x - 1) = T 1/2 (j) h (TX - r), 

^l(x) = r-^Mrx) - Mrx - 1). (33) 

Note that in the present case the following conditions are equivalent: 

A n G r _1 Z+ and A n+1 £ Z+ & \ n e r^A* & \ n e rZ+ = A+ (34) 

where we have denoted by A the generic elements of t _1 Z+. Consequently, 
the orthonormal basis of is given by: 

{<p s (x - \)} Xe(T 2 Z + +1) U {<p L (x - \)} XerZ + U {i) h (x - A)} AGrZ + , (35) 

and the orthonormal ("tau") Haar basis of L 2 (IR + ) is the set 

{r^L(r^-A)}. eZA6Tzr (36) 

The construction of the Haar wavelet basis corresponding to the substi- 
tution of the second type is carried out in a similar way. Lengths of tiles are 
respectively |L| = 1 and |S| = 1 — 1/(3. For a > 3 we have for the two scaling 
functions: 

<f>h(x) = l[o,i)(aO, (37) 
M% - a + 1) = (1 - l/P)-^ 2 l [a _ lia _ m (x). (38) 
Consequently, the refinement equations now read: 

a- 2 

Mx) = Yl ^ x - + (1 - l//9) 1/2 0s(/9x - (a - 1)), (39) 

1=0 

a— 3 

<k(x-a + l) = (l- 1/P)~ 1/2 MP* ~ Pifl - 1) - /) + 

1=0 

+ fa(/3x-(a-l)l3-a + 2). (40) 
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The corresponding orthonormal wavelet basis is constructed from the 
mother wavelet sets {ipL,i(x)}i=o and {ips,i( x )}iZo ■ Following the same method 
as in the previous case, we first pick the functions {(ft^x), 4>l(Px), 4>-l(I3x — 
1), . . . , (p^fix— (a— 1))} living on the interval [0, 1] and {(f>s(x), ^(flx), (fi^fix— 
1), . . . , — (a — 2))} living on [0, 1 — 1//3]. We then proceed to the Gram- 

Schmidt orthogonalization and L 2 -normalization. The orthonormal basis of 
Vi~ is obtained by collecting together all translates of the first orthonormal 
set at points of and all translates of the second orthonormal set at points 
of Ag . After removing functions in V + we get the basis of Wq . 

For the negative part of we just have to mirror the supports of scaling 
functions and wavelets. 

5.2 Non-orthogonal basis of Haar type 

We still have same scaling functions as in the orthogonal case. This means 
that all spaces are identical to the previous ones, but the basis of the 
orthogonal complement Wq of in is just required to be normalized. 
But the case j3 = r which leads to the same result as in the above, there are 
generically two mother wavelets ^ll and -^ls which live on tiles /9 _1 LL, resp. 
/3 _1 LS. The orthogonal complement of V$ in reads: 

W+ = span{^LL(:r - A LL ),^ LS (a: - A L s)} ALLe/3 -i A + L , ALse/ 3-i A + s , 

where Al L (resp. A^g) is the set of left-hand ends of words LL (resp. LS) in 

In the substitutional case (16) the refinement equations for wavelets are: 



Z+ 



^ll(x) = ' (MM - MP* - 1)), 

^ LS (x - (a - I)/ 0) = (j^lj ' (MP* " (a - 1)) " P 1/2 Mf3x ~ a))- 

The construction of the Haar wavelet basis corresponding to the substi- 
tution of the second type is carried out in a similar way. The corresponding 
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wavelets are given by: 



- (a - 2)/(5) 



f (a-l)P-l \ 

v 2/3-1 ; 




^sipx-a + l)). 



For the negative part of scaling functions and their equations are still the 
same but instead of the wavelet -?/>ls we now have to deal with a wavelet of 
the type Vsl- 



6 Lexicographical analysis of the Fibonacci 



We now focus on the rescaled version A = JF r /r 2 of the Fibonacci chain 
T T . This is a standard example of model set, well known by quasicrystal- 
lographers, and most of the obtained results about it are easily extendable 
to other cases described in Section 4. With the notation (22) this model set 
reads A = £[°' r ). We also view A = {A n } neZ as the increasing sequence of 
points X n < \ n+ i with Ao = 0. Let us give a list of important properties 
stemming from the interval nature of the window: 

• The model set A is self-similar with factor r 2 : r 2 A C A. More precisely, 
it is invariant under the amne-linear actions: 



• For a given n, A reads as a partition into subsets of n-letter words 
A ai a 2 ...a„ = {A | A is the left-hand end of the n-letter word a\a<i . . . a n , 
di E {L, S}} . More precisely, we have the following tree-like Fibonacci 
hierarchy of partitions: 



chain 



r 2 A = E [o,i) c A) r 2 A + S [0 ' r] = A. 



(41) 



A = 



A L U 
AiTuA^ U 



As, 



A 



A SL 



A = A LL s U A LS l U A SLL U A sls 
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x' 



x' 



Figure 4: Conjugate right nearest neighbour function / (left graph) and its 
iterated / o / (right graph). 

There are only n + 1 words built up from n letters (the substitution is 
Sturmian with minimal complexity) . 

• The above sequence of partitions is in one-to-one correspondence with 
partitions of the interval window [0,r 2 ). Accordingly, each element in 
the partition sequence is also a model set: A aia2 ... an = E^ 1,a ' 2 \ where 
u 1: u 2 e Z[t]. 

The properties of sets A aia2 ... 0ji , are suitably encoded by the one-to-one 
map / : A' — > A' introduced in [21] and defined by f(X' n ) = Ki+i where 
A n G A. This function uniquely determines the nearest right neighbour of 
any point A G A. Its reciprocal / _1 (A^) = \' n _i corresponds to the nearest 
left neighbour. The function / is explicitly given by: 



The graphs of the function f{x') and f(f(x')) are shown in Fig. 4. 

The function / divides the interval [0, r 2 ) into two parts. The first one 
is the interval [0, r), window of A L , and the second is [t, r 2 ), window of 
As- As one can see in Fig. 4 the next right neighbour of the conjugate r' = 
— 1/r of the discontinuity point r is the origin 0' = 0. Next, let us plot 
f^(x) = f(f(x)). This function encodes the distance to the second next 
right neighbour. The interval [0, r 2 ) is again split into two parts. This also 
means that 2-letter words have only two possible lengths. Generally, for each 
/N we will find a new point of discontinuity which splits the interval [0,r 2 ) 




x' + 1 if x G Al (since A n+ i = \ n + 1), 
x' — r if x G As (since A n+ i = A n + 1/r). 



(42) 
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into two parts and the n th next right neighbour of its Galois conjugate is the 
point 0. All points lying on the segment common to all iterates of / up to 
the n th one are the conjugates of left-hand ends of the same n-letter word 
- they induce a model set for this n-letter word. The following proposition 
easily ensues from this observation: 

Proposition 6. Left-hand ends of all possible n-letter words lie as the closest 
n left neighbours of the origin together with the latter, namely the points 
X-n, ■ ■ ■ , Ao- 

Proof. We know from the above that all these points are in different model 
sets and since there are n + 1 possible words we know that no one is missing. 

□ 

If we rearrange from the largest to the smallest one all n-letter words 
in the lexicographical order determined by L > S, A^ > > • ■ • > 

Ai n i > Ai n) , we find that A^ = E^'^i+V) for k = 0, 1, . . . , n, with A^ = 

£[4* n) >). Trie Aj^'s G {A_ n ,...,A } are left-hand ends of A^ n) and are 

ordered as A^ < Aj^ . 

Let us consider a n-tuple (A^, Afc+i, . . . , Afc +n ) G A n . We know that \k+i — 
Afe+i-i = 1 or 1/r. Consequently we have 



T 2 (\k+i - Afc + i_i - 1/r) = 



if \k+i and \k+i-i are ends of tile S, 

1 if \k+i and Afc+i-i are ends of tile L, 



which implies the following property. 

Property 7. The number ofL's in the word [A&, \k+ n ] is equal to 

n 

5^^ 2 (Afe+i - Xk+i-i ~ 1/r) = r 2 (A fc+ „ - X k - n/r). 
i=i 

Let us now characterize all points of £' 0,1 ) = r 2 ^' ' 1 " 2 ) in another way. 
From 



S^ 1 ) = 1 1 + kr \l, k G Z, I - - G [0, 1)1 =}► - < I < - + 1 =}► I = 
[ r J r r 

we easily recover a standard definition for chains of the Fibonacci type [22]: 
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Property 8. 



E [o,r») = _L E [o,i) = A 



A* 



+ kr) , k eZ 



It follows that the number of L's in any n-letter word is 



k + n 



+ (k + n)r 



kr — nr 



k + n 




~k~ 


T 




T 



k + n 


~k 






n 


k 


'k 






~n 




n 












+ -- 










or 




T 


T 






T 


r 


T 






T 




-T- 



Accordingly the number of S's in any n-letter word is n — |~^] = |_^J or 
n — [-J = ["-%] . In summary: 

Property 9. In any n-letter word there are |~^] L's and |_^J S's or |_^J L's 
and S 's. 



7 Fibonacci B-splines and wavelets 

7.1 B-splines scaling functions for the Fibonacci chain 

Using Def. 1 we have for s > 2, 

V^ S) (A) = {/ G C*- 2 , / G L 2 (M) | /| [An;An+l] is polynomial of degree < a - l}, 
where we remind that A n designates the n th element of A with An = 0. The 

( \ 

Riesz basis of Vq (A) is made up of translates of s + 1 functions (p^ , where 
the fij's, i = — s, . . . , 0, are all admissible s-letter words. Let us suppose that 
the support corresponding to the word is the interval [Aj, \+ s ]. Then we 
construct all s functions (p^ = <fii through the conditions below. 

• <j>i(x - Xi) = if x £ (Xi, X i+S ). 

• <f>i(x — Xi) is polynomial of degree s — 1 on the intervals [\i+k, Aj+fc+i], 
for < k < s - 1. 

• The (f>f\x — Aj)'s, < / < s — 2, are continuous in the points Xi + k for 
< k < s. 

• [ 4>i(x) dx = — . 

Jr s 
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7.2 Spline wavelets for the Fibonacci chain 

The construction of the wavelet basis of w'fl (A) rests upon the construction 
of the basis j^n' 5 '*, n E i?| of the space W^*\A). The first question to 

answer is how many different wavelet functions we have to build up and what 
are their supports. From (13) and Fig. 4 we see that A E = {A n | n E E} = 
{A„ E A | A n+ i ^ r 2 A} = Al- Hence we know that all supports start with 
L. The answer about the lengths of supports is afforded by Property 1. We 
recall that the support of the function ^k, k E E, is the interval (A&, ^k+N k ) 
which includes exactly 2s — 1 points from A\r 2 A and is the smallest possible. 
Let us now observe a few simple facts. We want to find all possible minimal 
supports containing 2s — 1 point of A\r 2 A. We know that any n-letter word 
supports n — 1 points of the Fibonacci chain. Since r 2 A = Ag + 1/t, the 
number of points of t 2 A in the interval (word) (A&, \k+ n ) is the same as the 
number of S's in the word (A&, \k+n-i}- If l^^r] oi Vyr\ is the number of 
points of t 2 A in the interval (A&, Xk+ n ) then [^^J or J 1 ^] is the number of 
points of A\r 2 A in the same interval. We now prove the following: 

Proposition 7. Let p be the length of the shortest word supporting 2s — 1 
points of A\t 2 A. Then all other minimal words starting with L supporting 
2s — 1 points of A\t 2 A have their lengths equal to p or p + 1. 

Proof. In words of length p there are \^-] = 2s — 1 or |_^J = 2s — 2 
points of A\r 2 A. The former words are precisely those we wish to find while 
we have to enlarge the latter with one point more. A successor (a (p + 1)- 
letter word which has the same p first letters) to words starting with L of 
the second type has to support 2s — 1 points from A\r 2 A also. Indeed, if it 
had 2s — 2 points from A\r 2 A only, we could remove the first L (and this 
removes also one point from A\r 2 A). Then we would have a p-letter word 
supporting 2s — 3 points from A\r 2 A and this is not possible. □ 

As it has been shown in the above, all words of length / = max fcg £ N k 
(from Proposition 7 we also know that 1—1 = min^g^ Nj,) can be found among 
those having left-hand end in the interval [A_/, Ao]. From the definition of Nk 
one can infer that there are 2s — 1 wavelets which overlap each other in the 
interval [A , Ai]. These wavelets are certainly all different - they have different 
supports. But the first wavelet the support of which starts from A_^ or from 
A_ i+1 (if A„2 + i E t 2 A) ends at A (it cannot end at Ai because A E r 2 A and 
also not at A 2 because then the length is larger than I). Consequently, we 
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have to select this wavelet also. Hence the total number of different wavelets 
is 2 s and we can finally assert the following: 

Proposition 8. For a given s there are exactly s + 1 B-spline scaling func- 
tions and 2s different wavelets. The supports of the latter have length equal 
to p = \(2s - 2)r] + f or to p + 1. 



7.3 Second-order splines for the Fibonacci chain 

From Fig. 3 and also from the underlying substitution rules, we can see that 
for s — 2, there are just 3 possible words, namely LL, LS, SL. Hence we get 
the following scaling functions: 



x s e [o,i], 

2-x x e [1,2], 



f x + t xe[-r,-l/r], 
1 + ] \r 2 -r(x + r) xE [-l/r,0], 



<J) S l(x + 1/t) 
The scaling equations are: 



t(x + 1/t) i6[-l/r,0], 
r- (x + 1/t) x G [0,1]. 
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<Pll(x) = \<Pll{t 2 x) + ^<p LS {T 2 x - 1) + <Psl{t 2 x - 2) + 

-0ll(t 2 x - r 2 ) + ^0Ls(r 2 x - r 2 - 1), 
r r° 

0Ls(a: + r) = ^0 LL (r 2 a: + r 3 ) + ^</) L s(t 2 x + r 2 + 1/r) + 

0sl(t- 2 ^ + r + 1/r) + -^0ls(t 2 x + r), 

SL (x + 1/r) = -(p LS (r 2 x + r) + SL (r 2 a: + 1/r) + -<p hh ( T 2 x) + 
r r 

l0 LS (r 2 x-l). 
r° 

So in this case the Riesz basis of (A) is the set of functions 

{0ll (x - All), 0ls(^ - A LS ), 0sl(^ - AsL)} ALLe A LL ,A LS eA LS ,A SL eA SL ■ 

In order to make the explicit computation of wavelets more transparent 
we recall in Table 1 below the points of the Fibonacci chain A = {X n } n& z 
around zero. 



L 


L 


S 


L 


S 


L 


L 


S 


L 


L 


S 


A-5 


A-4 


A-3 


A-2 


A_i 


Ao 


Ai 


A 2 


A 3 


A 4 


A 5 


-r 3 


-r 2 -± 

T 


-T-i 

T 


— r 


_ i_ 

T 





i 




r 2 


r 2 + l 




-1-2t 


-2t 


l-2r 


— T 


1-T 





i 


2 


1 + T 


2+r 


3+r 


yes 


no 


no 


yes 


no 


yes 


no 


no 


yes 


no 


no 



Table 1: Fibonacci chain points located around the origin. 



The first row tells us whether the point is a left-hand end of the tile L or S, 
the second row is just the indexation of points, in the third row there points 
of the Fibonacci chain are written in terms of r-expansion, in the fourth row 
they are written in the form a + br. The last row indicates whether the point 
is element of r 2 A or not. 

Let us now construct the wavelet basis of Wq 2 \A) by applying the re- 
sults of Section 3. We look at first for the ^n S ~^ 's, n G E, which form 
the basis of W^i(A) = |/ G | /| r 2 A = oj . Their respective supports 

are [A„, A n+ jvJ, where n G E = {n G Z | A n+ i ^ r 2 A} and N n is the 
smallest number for which the equality N n = 4 + #{(A n , \ n+ N n ) H r 2 A} 
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Cllsls(2 


-A_ 5 ) 


Clslsll(^ — A_ 4 ) 


tile 


interval 




Qi 


h 


Hi 


L 
L 
S 
L 
S 
L 
L 


\X X d) 

Y 05 /x — 4/ 

[A_4,A_ 3 ) 

[A- 3 ,A_ 2 ) 

[A_2,A_i) 

[A-i,A ) 

[Ao, Ai) 

[Ai,A 2 ) 


g 

6(l+26r) 
11 

12(18+17r) 
11 

24(l+4r) 
11 

12(3+r) 
11 


o 

6 

12(5-13r) 
11 

12(4+5r) 
11 

12(2-3r) 
11 


6 

-24(l+2r) 

6(91+123r) 
11 

6(127+224r) 
11 

21(13+16r) 
11 

3(13+16r) 
11 



6 

-6(3+4t) 

6(58+79t) 
11 

18(13+16t) 
11 

3(13+16t) 
11 












norm 


£ ^6(270+431t) =13.8519 


£ v /3(16754+27145r)=77.572 



Table 2: Analytical expression of Cllsls(^ - A_ 5 ) and Clslsll(z - A_ 4 ). 



holds true. We recall that the equality A E = {\ n £ A \ n <E E} = A L 
determines that all supports will start by L. The set of these supports is 
= {LLSLS, LSLSLL, LSLLS, LLSLL}. This means that we have four 
different shapes of s, one per support. As second derivatives of func- 
tions ^n\x), n G E we obtain four different functions (\ n {x), precisely 
( Xn (x - A„) = j^Wn\x), n G E (see Fig. 6). Due to the equality (\ n (x) = 
(\ k (x) if A n , Afc G A M where ji G Q^, we just denote the latter by (^x). 
From Table 2 and 3 we deduce the analytical expressions of functions 




Figure 6: Four cubic functions (x)'s and their second derivatives 
Cx n (x - A n )'s for n G {-5, -4, -2, 0}. 
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Clslls(^ 


-A_ 2 ) 


Cllsll(^ — A ) 


tile 


interval 


ki 


% 


h 




T 

L 

s 

L 
L 
S 
L 
L 


[A-2, A_iJ 

[A-i,A ) 

[A ,AO 

[Ai,A 2 ) 

[A 2 ,A 3 ) 

[A 3 , A 4 ) 

[A 4 , A 5 ) 


6 

-24(1+2t) 
3(14+19t) 
-3(10+19t) 
6(4+5r) 




6 

-6(3+4t) 
3(8+11t) 
-6(1+4t) 


6 

6(-25+42r) 
11 

36(10+3r) 
11 

42(-1+3t) 
11 

6(-l+3r) 
11 




6 

36(-6+7t) 
11 

36(-1+3t) 
11 

6(-1+3t) 
11 












norm 


^/6(191+309r)=64.3882 


^ v /6(281+411t)=13.6981 



Table 3: Analytical expression for Clslls^ — A- 2 ) and CllsliX^ — A ). 



jj, G Vt^, translated to their basic admissible intervals. On each subinterval 
[Aj, Aj + i) the function A M ) is described by a linear function ki(x— Aj)+^. 

The values of fcj's and g^'s are given in these tables. 

Remark 2. It should be noticed from all elements of the extension field Q(r) 
shown in the tables that imposing \I/^(A n+ i) = 1 is not an optimal choice. 
The simplest expressions are found with the choice \l4 4 ^(A n+ i) = 11/3. The 
appearing of number 11 is due to an interesting algebraic feature. In order 
to determine the coefficients of the linear ( 's we have to solve a system of 
equations with coefficients in the field Q(r), and so all solutions are of the 
form l = ^g = {a+ ^- + d J d T) E Q(t). Now, it happens that the "algebraic 
squared norms" ww' = c 2 — d 2 + cd of denominators are equal to ±1 (i.e. w 
is unit in 7i[t\) or to ±11. More generally, properties and shapes of wavelets 
are independent of the fixing ^n\\ n+ i) = u. With the generic latter choice, 
all values in the tables 2 and 3 should be multiplied by u. 

By dilation and normalization of the (^(x)'s we eventually obtain four 
mother wavelets ip^(x) = ^fe-^i , see Fig. 7 for their basic admissi- 

ble translations. Then the set of functions — ^en i ,,K l _ l er- 2 A IJ , forms 

(2) 

a Riesz basis of Wq (A) and consequently the set of functions 
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forms a Riesz basis of L 2 (R). 




7.4 Construction of scaling equations for wavelets 

Let us here describe the algorithm for obtaining scaling equations in the case 
of piecewise linear wavelets. For convenience we first set up the equations for 
the functions \i G O^. The corresponding equations for wavelets ip^, jj, G 
Qip are trivially derived from them. Let us denote by Q^, the set of words 
supporting the scaling functions (in the present case = {LL,LS, SL}). 

(2) (2) 

Due to W_{ C Vq we can write 

CO - V) = 9u,iM x ~ 0- 

The translate to the point Aj G A u of a scaling function (p v for all G 
reads as: 

£ - A .; 



Ai (x - Ai) 




for x G [Ai, A i+ i] 

Aj+2 — x 



A,:+o — Ai_n A,:_lo — Aj 



for x G [A i+ i, A i+2 ]- 



We look for g^'s in the scaling equation 

i+Ni-2 

( Xi (x-K)= Yl 9i^-\i). (43) 

j=i 
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Because of support condition and continuity at the point Aj we can write 

ki(x - Xi) + qi\ x=x . = 0, 

and so the coefficient qi has to be 0. Hence we have on this interval the 
equation 

ki(x - Xi) = g~ *— =>• gi = ki(X i+1 - Xi). 

The function is continuous at Aj+i so we have 

h{x - \)\ x=Xi+1 = k i+1 (x - A i+1 ) + q i+ i\ x=Xi+1 , 

thus 

q i+ i = h(X i+1 - Xi) = (44) 

• [Ai+i, Aj_|_2] : 

On the interval [X i+ i, X i+2 ] we have the equation 

/ \ A^_|_2 — ^ ^ — A^i 
k i+1 (x - X i+1 ) + q i+1 = g- h g i+ i- . 

Using the equation (44) we obtain g i+1 = k i+ i(X i+ 2 — X i+ i) + q i+ ±. From the 
continuity at the point A« + 2 we also have 

h+i(x - Aj+i) + qi+i\ x=x . +2 = k i+2 (x - X i+2 ) + 9i+2| x=A . +2 , 



thus q i+2 = g i+ i. 
• [Ai+j, Ai_|_j_|_i] : 

By induction we have generally for /, < I < Ni — 2 

g i+ i = k i+ i(X i+ i + i - X i+t ) + q i+ i = q i+ i+i. (45) 

This method is simple and very general. Actually it can be carried out for any 
self-similar locally finite Delone set in the case of piecewise linear wavelets. 

Finally the scaling equations for basic admissible translations of functions 
i/jn, \i G have the form 
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^LLSLs(^-K-s) 


V>LLSLs(z-K-4) 


</>lslsll(z-K-2) 


^llsll(^) 


<Pll{t 2 x - A_ 5 ) 


6 








(phs(r 2 x - A_ 4 ) 


g(5-13r) 


6 






<Psl(t 2 x - A_ 3 ) 


i(4 + 5r) 


-6(3 + 4r) 






<hs(r 2 x - A_ 2 ) 


g(2-3r) 


g(58 + 79r) 


6 




<Psl(t 2 x - A_i) 




-g(13 + 16r) 


-6(3 + 4r) 




<Pll(t 2 x - A ) 




g(13 + 16r) 


3(49 + 79r) 


6 


0Ls(r 2 a; - Ai) 






-6(1 +4r) 


g(6-7r) 


<Psl(t 2 x - A 2 ) 








-g(l-3r) 


0LL(r 2 x - A 3 ) 








g(l-3r) 



Table 4: Table of coefficients for scaling equations of ip^, ji G f^. 



Coefficients g^ t are recalled in Table 4. The norms in denominators are: 

IICLLSLs(r 2 x)|| 

||CLSLSLL(r 2 a;)|| 
IICLSLLs(r 2 x)|| 
IICLLSLs(r 2 x)|| 

8 Wavelets for stone-inflation tilings 

The construction of spline wavelets for self-similar Delaunay point sets in 
M. d has also been considered by Bernuau in [8]. Here we shall briefly de- 
scribe similar constructions for a tiling T of M, d which has the stone-inflation 
symmetry aT C T where a = (b,9R) E R d x (R* + x SO(d)), 9 > 1, is an 
affine transformation of M. d . Related constructions were also described in 
[23, 24]. Let us first recall what we understand by tiling. A tiling of M, d is a 
covering with non-overlapping pieces all congruent (in the Euclidean group 
action sense) to tiles belonging to a predetermined finite set of prototiles or 



= — ^6(109 + 161r) = 8.561, 

= ^^3(6363 + 10391r) = 47.942, 
= ^6(73 + 118-7-) = 39.794, 
= -^v 76 ! 151 + 130r) = 8.466. 
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Figure 8: First three steps of stone-inflation for a Penrose tiling. 



tile alphabet, A = \T\, . . . , T n }: 

K d =|J P *' p i C °- T ^ T i^ A > Pi^PjCR d ', t^j,d'<d. (46) 

-ieN 

In this definition, a tile is supposed to be compact, equal to the closure of 
its interior, and homeomorphic to a topological ball. Moreover, in the tiling 
any two tiles have pairwise disjoint interiors (all tiles of the tiling pack IR d ). 
Many tilings T of M. d have the so-called stone-inflation symmetry: all tiles 
of T, when rotated and scaled by a given OR > 1, 9 > 1, R E SO(d), and 
translated by a given b G M. d , can be packed face-to-face from the original 
ones. An example is provided by a tiling in the Penrose class as shown 
in Fig. 8. Another example, also relevant to quasicrystalline studies is the 
three-dimensional Danzer tiling [25]. 

More precisely, suppose we are given a tiling T of W d built from a finite 
set of prototiles {T 1; . . . , T n } present in the tiling, 

n 

T = U P * = U U Ti- T " ( 4? ) 

ieN ;=i 7i gr ; 

where 7^ = (fy, Ri) is a translation-rotation in the Euclidean group R d xiSO(d) 
and C M d x SO(d) is made of all those transformations 7/ = (bi,Ri) 
(including the identity) which bring the prototile TJ to one of its congruent 
companions appearing in the tiling. Here, we keep the freedom to restrict 
the r^'s to pure translation sets, at the price of enlarging the number of 
prototiles. The stone-inflation symmetry based on the affine-linear inflation 
a = (b, OR) then means that for each / and each 7/ 6 T h the following finite 
patch 

n 

aPi = a lx ■ T x = |J |J a lm ■ T m , (48) 
m=i o- Jm er m 
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is present in the tiling. Hence, it becomes possible to deal with an infinite 
sequence of inflated-deflated versions of the tiling T 



• • • C a~ j T = [J <j~ j ■ P t = T 3 (ZT j+1 ..-. (49) 

ieN 

Here, the inclusion relation between two tilings should be understood as: 
T C T' if any tile of T is patch of tiles of T'. Let A(T) denote the set of 
vertices of the tiling T. It is clearly Delaunay and we have the counterpart 
of (49): 

•••A^cA^cA^).... (50) 
Furthermore, we have the denseness property: 

R d = \jMTj). (51) 

j 

8.1 Haar wavelet basis 

In order to prepare the discussion about the feasibility of finding spline 
wavelets adapted to tilings of the above type, let us first adapt to our context 
the construction of an orthonormal Haar wavelet basis of L 2 (M. d ) proposed 
in [8]. We first define a sequence of spaces (Vj(T))j e z such that each space 
Vj(T) is the closed subspace of L 2 (M, d ) of functions which are constant on all 
tiles <?~ j Pi, i6N. Then we have, 

Proposition 9. The sequence (Vj(T))j e z is a cr-multiresolution analysis of 
L 2 (M. d ) in the following sense: 

(i) for oHjeZiVj-WcVjiT), 

(ii) {J jeZ Vj(T) is dense in L 2 (R d ), 

(m) n je z^cn = {o}, 

(iv) for all j e Z,x G R d , /(x) G V 3 {T) ^ f(a~^) G V (T), 

(v) there exists a finite number n of scaling functions, 0i(x), 2 ( x ), • • • , 0n( x ) 
such that all their admissible linear-affine transformations {(fiii^ 1 • 
x)}i<Kn,7 ; er; form an orthonormal basis in Vq{T). 

Proof. 
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(i) This inclusion results from the stone-inflation property of tiling T, 
aTdT. 

(ii) This is true through the fact that every continuous function / with 
compact support on M d can be written as uniform limit of the sequence 
(fj)j>o, 

/i( x ) = Yl /( X i>™) 1 a- J p m (x), 

where x j>m G o~ J P m . 

(iii) By construction it is clear that only the function /(x) = is included 
in all spaces Vj(T). 

(iv) Let us choose /(x) G Vj(T). Then we have 

/( x ) = c m l CT - 3Pm (x), with \ c m\ 2 < °°> 

rreSN m£N 

and this is equivalent to 

/(<7 -J x) = ^C m l Pm (x). 

Thus /(<j-'x) G Vb(T). 

(v) We have n scaling functions 0i(x), fa (x), . . . , 0„(x) which are the nor- 
malized characteristic functions of the corresponding tiles TJ, / = 1, . . . , n 
divided by the square root of its volume, e.g. for a tile X <E A 

X (X) = l x (x)/y/\X\ 

where \X\ means the volume of the tile X. Due to normalization we see 
that all admissible linear-afline transformations of n functions {</>;}™ =1 
form an orthonormal basis of Vq(T). 

□ 

We now proceed to the construction of the Haar wavelets. For each 
prototile 7], 1 < I < n there are finitely many tiles P^, such that: 

aTi = P h U • • • U P lki -<=>• Ti = a- x P h U • • • U a^P^ , 
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and the set |P/ 1; . . . , p fc; j forms a patch present in the tiling T. Denote 

by Vo^T) the subspace of V (T) of functions which are zero (almost every- 
where) outside the tile Pj. Accordingly we define Vj tPi (T) as the subspace 
of Vj(T) of functions equal to zero (almost everywhere) outside the tile Pj. 
Therefore we obtain the orthogonal decompositions: 



V (T)=($ ± V , Pi (n (52) 

VjCn=®±v s An (53) 



and the inclusions 

K,,p t (T) C V 1:Pi (T) c • • • C lfo(T) • • • . (54) 
The wavelet space Wo(T) is the orthogonal complement of Vq{T) in Vi(T), 

Vi(T) = V (T) © ± W (T). (55) 

More generally 

We also define the orthogonal complement of Vj :Pi (T) in Vj +1:Pi (T) as 

^ +1 , Pi (T) = V hPi {T) ©j_ W jiPi (T). (56) 

It is clear that 

^■co = 0±wwn 

Finally we get the orthogonal decomposition of the whole Hilbert space: 

L\M. d ) = 0j_WO(T). 

Thus the construction of Wo(T) is equivalent to the construction of all 
Wo tPi (T)'s. Since any tile Pj can be written as 

P = 7^ = Jiia-'Pt, U • • • U (T -1 ^), 

it is sufficient to find wavelets for prototiles 7], I — 1, . . . , n, and the whole 
basis of W / (^') w ih be formed by all admissible linear-afline transformations 
of these "protowavelets" . There results the following proposition. 
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Proposition 10. For every prototile T h I — 1, ... ,n, given also by 
we have ki — 1 orthonormal wavelets 

Proof. Let us denote by Vo,^ (X) the space of functions constant on tile TJ and 
equal to zero otherwise, i.e. those functions proportional to 1t ; (x). We then 
denote by Vij^T) the space of functions constant on tiles a~ x P\^ . . . , a~ l Pi ki 
and otherwise equal to zero. The space of wavelets corresponding to the tile 
Ti is found as the orthogonal complement of V^^iT) in Vi,t ; CO, i-e. 

Vm{T) = V m (T)®±W m {T). 

One possible construction of an orthonormal basis of Wo,T t {T) is as follows: 

we start from the set of functions {1t ; (x), 1 ct -i p ; (x), . . . , Pl (x)} which 

i fy— i 

form a basis of Vi^T). Because of the presence of the function lr,( x ) we 

have put aside (for instance) the last function l a - 1 P l ( x ). We then proceed 

k i 

to Gram-Schmidt orthogonalization and L 2 -normalization. Thus we get an 
orthonormal basis of Vi^T) and after removing the single basis element in 
Vq ^ {T) we are left with ki — 1 functions which form an orthonormal basis of 

wwn □ 

The orthonormal Haar basis of L 2 (M. d ) adapted to the tiling T is finally 
given by: 

U U U { 93i ^ (7f V • x) , . . . , ^1^, ( 7l - V • x) } . 
jezi<i<n 7; er( 

8.2 Beyond Haar: spline wavelet basis 

Let T be a stone inflation tiling of M, d . Suppose prototiles are polytopes and 
the set A(T) of vertices is of finite local complexity. How to build piecewise 
linear, compactly supported, spline scaling functions and related wavelets 
living on tiles of T? More precisely, one could think about going through 
the following steps. 

• How many tiles are needed for building the minimal patches supporting 
scaling functions? 
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• How many different patches of such type exist in the tiling? 

• How to characterize such spline functions from a functional analysis 
point of view? 

• Suppose the previous questions answered. How to build the correspond- 
ing multiresolution? 

• Suppose the previous question be solved and related scaling functions 
be determined. How to find the corresponding wavelets? 

• Are these wavelets compactly supported? 

• If yes, on which patches in the tiling? 

• How to extend this material to smoother splines? 

It is not the aim of the present paper to achieve this ambitious pro- 
gram. Let us just sketch which procedure could be followed in the simplest 
two-dimensional case. In order to build "pyramidal" functions (the two- 
dimensional analogue of the one dimensional spline "hat" functions) which 
are to play the role of scaling functions, we have to accomplish a triangulation 
of the tiling T (like we already have for Penrose tilings). By triangulation we 
mean that in each prototile some of the vertices could be pairwise connected 
with a new segment (in order to divide the prototile into triangles) under 
the condition not to create new nodes (the new edges should not cross). All 
tiles congruent to the divided prototile will be triangulated in the same man- 
ner. Then to every point A G A(T) there is an associated pyramid. The 
support of this pyramid is delimited by the neighbours of A. (A neighbour 
is a point connected with A by an edge). The vertex of the pyramid is obvi- 
ously located right above A. Because of the finite local complexity we know 
that there exists a finite number of different pyramids. These pyramidal 
functions are uniquely determined once a certain normalization condition is 
fulfilled. These pyramids and all their admissible linear-afline transformation 
then define the space Vq(T). The multiresolution analysis follows. 

A simple example is provided by the square lattice A = Z 2 : we just have 
to divide each square in the same manner into two triangles. Then we get 
one scaling function which is a "hexagonal pyramid", as is shown in Fig. 9. 

Another example is provided by a five-fold Penrose tiling. Seven pyra- 
midal functions exist here. As is shown in Fig. 10, their respective supports 
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1 2 3 4 5 6 



Figure 9: Support of "hexagonal" pyramidal function for the square lattice 
A = Z 2 . 

are (in terms of kites, darts and rhombuses): a decagonal patch of 5 kites, 
a star-shaped patch of 5 darts, a kite-shaped patch of 2 kites and 1 "fat" 
rhombus (i.e 2 obtuse triangles), a rhombus-shaped patch of 1 dart and 1 
kite (an "ace"), a patch of 3 darts and 1 "thin" rhombus (i.e 2 acute trian- 
gles), a patch of 1 dart and 2 thin rhombuses, and a patch of 3 kites and 2 
obtuse triangles. 

9 Conclusion 

The Bernuau construction of wavelets adapted to self-similar aperiodic point 
sets has been carried out in dealing with some one- dimensional examples. 
Our aim is currently to apply these wavelets to the analysis of aperiodic 
structures, like diffraction spectra of Fibonacci chain, and to compare our 
results with more standard wavelet analysis (e.g. dyadic wavelets). Of course, 
an essential step in decomposition and recomposition schemes will be the 
determination of corresponding biorthogonal basis. We also plan to extend 
our constructions to higher-dimensional cases, like Penrose or Danzer tilings, 
in view of practical applications to quasicrystalline studies. 
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Figure 10: Supports of seven different pyramidal functions for a (five-fold 
symmetrical) Penrose tiling. 
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